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Abstract 

The dynamics of a glass-forming material slow greatly near the glass transition, and 
molecular motion becomes inhibited. We use confocal microscopy to investigate the 
motion of colloidal particles near the colloidal glass transition. As the concentration 
in a dense colloidal suspension is increased, particles become confined in transient 
cages formed by their neighbors. This prevents them from diffusing freely through- 
out the sample. We quantify the properties of these cages by measuring temporal 
anticorrelations of the particles' displacements. The local cage properties are related 
to the subdiffusive rise of the mean square displacement: over a broad range of time 
scales, the mean square displacement grows slower than linearly in time. 
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1 Introduction 



As glass-forming materials are cooled, the sharply increasing viscosity of the 
liquid is accompanied by equally dramatic changes in the motion of tracer 
particles within the material [1]. In particular, the mean square displacement 
(Ax 2 ) (MSD) of an ensemble of tracer particles embedded in a glass-forming 
material forms a plateau at intermediate lag times, reflecting the crowding of 
the particles which prevents easy rearrangements [see Fig. 1(a)]. At longer lag 
times, the MSD shows an upturn, returning to diffusive motion, albeit with 
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a greatly reduced diffusion coefficient ((Ax 2 ) ~ At). Dense colloidal sus- 
pensions are simple materials which undergo a glass transition as the particle 
concentration increases, and provide a way to directly study the anomalous 
kinetics of the colloidal particles near the glass transition, to determine how 
the local motion of individual particles gives rise to the unusual behavior of 
the ensemble MSD [2,3]. The plateau in the MSD is subdiffusive: for a range of 
time scales At, (Ax 2 ) grows as (Ax 2 ) ~ (At) 7 with 7 < 1; 7 = 1 is the more 
typical diffusive case. Typically sub diffusion arises when a system possesses 
memory [4]. In this work, we test this by looking for temporal correlations 
in the particle motions. We find that these correlations do exist and are due 
to the "cage effect" of glassy systems (see Fig. 2). We characterize this cage 
effect, and directly connect the local description of particle caging to the sub- 
diffusive plateau in the MSD and the lag-time dependent anomalous diffusion 
exponent 7(At). 



2 Experimental procedure 



Our samples are colloidal poly-(methylmethacrylate) (PMMA) particles, ster- 
ically stabilized by a thin layer of poly-12-hydroxystearic acid [2,5,6]. They 
are in an organic solvent mixture of cyclohexylbromide and decalin, chosen 
to closely match the density and index of refraction of the particles [5]. The 
particles have a radius a = 1.18 /im and a polydispersity of ~5%. They are 
dyed with rhodamine dye, which results in a slight charging of the particles. 
Despite this slight charge, their phase behavior is similar to colloidal hard 
spheres [7]: we find 0f reez e = 0.38 and me i t = 0.42 (for hard spheres these val- 
ues are 0/ = 0.494 and <p m = 0.545). As the concentration is further increased, 
we see a glass transition at <j) g ~ 0.58, in agreement with what is seen for hard 
spheres. Samples with > 4> g do not form crystals in the bulk even after they 
have been sitting for several months. Moreover, the diffusion constant for such 
samples goes to zero - the samples become nonergodic [1]. 

We view the colloidal particles with a fast scanning laser confocal microscope, 
to obtain three-dimensional images from deep within the sample [2,3,5,8]. In 
practice, we focus at least 30 /xm from the coverslip of the sample chamber, to 
avoid wall effects. By taking a series of three-dimensional images at intervals 
of 10-20 s, we are able to follow the motion of several thousand colloidal 
particles for several hours. We identify particle centers with an accuracy of 
0.03 /im horizontally and 0.05 /im vertically; the poorer vertical resolution is 
due to optical limitations of the microscope. For further details, see Refs. [5,9]. 
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3 Results 

We calculate the mean square displacement (Ax 2 ) from the measured particle 
positions, and several typical curves are shown in Fig. 1(a). The data at short 
At (less than 10 s) is obtained from two-dimensional measurements within 
the three dimensional sample, in order to improve the time resolution. At the 
shortest lag times, (Ax 2 ) increases due to the diffusive motion of the particles. 
At intermediate time scales, the MSD has a plateau, which becomes more 
pronounced as the volume fraction <f> increases toward <fi g ~ 0.58. This plateau 
is due to the cage effect: particles are trapped in transient cages formed by 
their neighbors, and thus cannot diffuse freely through the sample [1,10]. At 
the largest At, the cages rearrange, and particles are able to diffuse throughout 
the sample, albeit with a greatly decreased diffusion coefficient [10-12]. 
This can be seen in the particle trajectories shown in Fig. 2. Figure 2(a) shows 
two-dimensional projections of trajectories of several particles within a small 
region. The two particles marked b and c are magnified to the right, and show 
the difference between caged motion, and the rearrangements. 

The cage rearrangements - that is, the relatively rapid shifts in particle po- 
sitions seen in Fig. 2 - are reflected in broad tails for the distribution of 
particle displacements [12-15]. These distributions are shown by the symbols 
in Fig. 3(a). The time scales for the displacements are chosen to be compa- 
rable to the end of the MSD plateau. The majority of particles move only 
short distances, as they are confined within cages. However, the distributions 
show that a nontrivial fraction of particles do move large distances, more than 
would be expected if the distributions of displacements were gaussian [dotted 
lines in Fig. 3(a)]. A traditional way to quantify the relative size of the tails 
of the distribution is to calculate a nongaussian parameter, which compares 
the fourth moment of the distribution to the second moment: 

, A , 3(Ar 4 (At)) 1 

which is zero for a gaussian distribution, and larger when the distribution is 
broader (for example, a 2 — 1 for an exponential distribution) [12,15-17]. This 
parameter is close to zero at small and large lag times At, and is a maximum 
at an intermediate value At* which we use to define the cage rearrangement 
time scale [11,12]. This time scale is indicated by vertical bars in Fig. 1(a), 
and corresponds qualitatively with the end of the MSD plateau. 

To quantify the cage effect, we wish to look for temporal correlations in a 
particle's motion; we follow the method of Doliwa and Heuer [10,11]. In par- 
ticular, if a particle moves in one direction for a period of time, its neighboring 
particles (the "cage") will prevent further motion in that direction, and may 
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Fig. 1. (a) Mean square displacement for three "supercooled fluids," with volume 
fractions <j) as indicated. The vertical lines indicate the cage rearrangement time 
scale At*, (b) The symbols indicate the measured anomalous diffusion exponent 
J (At), equivalent to the logarithmic slope of (Ax 2 ). The lines show the predicted 
value 7est based on Eq. 3 (dotted line <f> = 0.46, solid line 4> = 0.52, and dashed line 
(f> = 0.56). (c) The anticorrelation scale factor c(At) from Eq. 2; see text for details, 
(d) r cagc as a function of At. The symbols in (b-d) are the same as part (a). 

push the first particle back toward the middle of the cage. In this way, the 
positions of the neighboring particles, which have shifted slightly to allow the 
interior particle to move, provide a "memory" of the interior particle's motion. 
Thus we expect that usually a particle's motion will be temporally anticorre- 
lated, unless it is involved in a cage rearrangement (and thus moves and then 
stays in its new position). To look for this, we pick a time scale At and then 
consider displacement vectors for each particle, Ar mn = r(nAt) — r(mAt). 
In particular we wish to determine how Ar*i 2 depends on Ar i, and how this 
depends on the time scale At [10]. 

Anticipating that Ar*i 2 is directionally correlated with Afoi, we consider two 
components of Ar*i 2 : x 12 is the component of Ar 12 parallel to Ar i, the original 
displacement, and y i2 is the component of Ar 12 along an arbitrarily chosen 
direction perpendicular to r 01 [10]. Because of the arbitrariness in calculating 
yi2, the average (yu) = 0. For dilute samples, caging does not occur and 
(^12) — 0; (X12) will be negative if memory effects are present. Particles which 
initially move farther must move their neighbors farther as well, and so we 
expect that (x\ 2 ) will depend on how far a particle has originally moved, 
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Fig. 2. (a) Trajectories of particles from a sample with 4> = 0.52, over a 2 hour 
period. The axes are labeled in microns, and the circle illustrates the particle size. 
These trajectories are from particles within a 2.5 /im thick region within the sample; 
the gray shades indicate vertical distance (darker is closer to the coverslip). (b) and 
(c) are magnifications of two of the trajectories, with tick marks indicating 0.2 //m 
spacings. These two particles alternate between being trapped in a local cage, and 
a slight jump to a new location when the cage rearranges. 
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Fig. 3. (a) Probability distribution functions for displacements roi with time scales 
At = 260 s for 4> = 0.46, 700 s for cf> = 0.52, and 1000 s for <p = 0.56. (b) (x 12 ) as a 
function of roi for the same data shown in (a). The values of At for the three data 
sets have been chosen to produce similar behavior at small roi, which in these cases 
is reasonably well described as (xu) = — (0.26)roi (indicated by the dashed line). 
The departure from the small roi behavior occurs at r cage ~ 0.75, 0.35, 0.25 fxm. for 
= 0.46,0.52,0.56. 

r oi = \tqi |. To investigate this we compute the average value (£12) as a function 
of roi, and plot this in Fig. 3(b) for three different volume fractions [10]. At 
has been chosen so that the curves have similar behavior at small r i, and also 
to be close to the cage rearrangement time scale At*. The average is taken 
over all particles and all initial times, xu is negative, indicating anticorrelated 



5 



motion: particles which move in one direction during the first time interval 
will, on average, move in the opposite direction during the subsequent time 
interval. This is a direct signature of the cage effect. Moreover, for particles 
with small displacements r i, the average subsequent displacement (xi 2 ) is 
linearly proportional to r i, as indicated by the dashed line in Fig. 3(b). For 
larger r i, (^12) is no longer proportional to r i, and in fact becomes almost 
independent of r m [10]. 

The departure from the linear behavior at small r 01 occurs at smaller distances 
as the volume fraction increases toward the glass transition. The existence 
of two regimes - a linear response at small r m and a breakdown of this linear 
response at larger tqi - suggests that the crossover point can be taken as r cage , 
and the two regimes be identified as caged particles and rearranging particles 
respectively. In other words, particles with r 01 < r cagc typically remain caged, 
and the effect of the cage is to push the particle back toward its original 
position [10,11]. The strength of this effect is given by 

(^12) = -cr i, (2) 



with for example c = 0.26 for the data shown in Fig. 3(b). Particles with 
r oi > r ca g c still tend to be pushed back, but not as far as predicted from linear 
extrapolation from the small r i behavior: thus these particles may end up 
in new positions, and their behavior reflects cage rearrangements rather than 
caged motion. The changes seen in Fig. 3(b) as <fi is increased shows that the 
cage size r cage decreases as the glass transition is approached [18,19]. 

By studying the At dependence of the proportionality constant c and cage 
size r cage , we can better understand the MSD. The value of c depends strongly 
on the chosen time scale At, as shown in Fig. 1(c). In the middle of the MSD 
plateau, c is large, close to 0.5; at larger At it decreases, signaling a diminishing 
cage effect. c(At) can be related to the logarithmic slope of the MSD [10], to 
directly connect the cage effect to the subdiffusive MSD plateau. Locally the 
MSD grows as (Ar 2 ) ~ At 7 ^ A ^, with the anomalous diffusion exponent 7 (At) 
equal to the logarithmic derivative of (Ax 2 ). This can be estimated as: 



, rfln(Ar 2 ) 
Test (A*) = -JJ^f- 

_ ln[|Ar i + Ar l2 \ 2 /(r 2 m )] 

ln(2At/At) 
_ ln(2 + 2(:r 12 ro 1 )/(r 2 1 )) 
ln2 

«l + ln(l-c(A*))/lii2. (3) 
We have used (r 2 2 ) = (r^) (time invariance) and the final approximation uses 
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Fig. 4. ctii (connected symbols) and <tj_ as a function of roi, for three different volume 
fractions as indicated. The time scales are as in Fig. 3. 

(x 12 r 01 )/ (7-57) w (xi2) / (r i) ~ — c, in analogy with Eq. 2; we have verified that 
these approximations are reasonable [11]. In Fig. 1(b) the symbols show 7 (At) 
computed directly from the MSD, and the lines show 7 es t(At) calculated from 
Eq. 3. The subdiffusive plateau in the MSD is seen as a broad range of At 
for which 7 (At) < 1, although it is also clear that 7 does not have a constant 
value anywhere in the plateau, but rather is a smoothly evolving function of 
At. Moreover, the behavior of 7(At) is well-captured by the calculated value 
based on c(At), as shown by the agreement between the symbols (7) and the 
lines (7 C st)- In other words, the subdiffusive behavior of the MSD is a direct 
consequence of the caged motion of the particles, as measured by Eq. 2. As 
At increases, the cage effect becomes less important, c(At) decreases toward 
zero (no caging), and the MSD approaches diffusive behavior (7 — > 1). 

The behavior of the cage size r cagc is shown in Fig. 1(d). r cage is relatively 
insensitive to At, indicating the the size of the cage is more likely a static 
property [10,18]. The cage size decreases as the glass transition is approached, 
although it has a nonzero value at the glass transition [18,19]. The diffusive 
behavior of the MSD at large time scales can thus be thought of as due to 
the random walks of the individual particles, each taking steps of size r cage in 
random directions [18]. On the cage rearrangement time scale At*, only a few 
particles move (5 - 10 %) [2,12,18], and so in fact the average time between 
random walk steps is much larger than At* as seen in Ref. [18]. 

Further insight into the cage effect can be found by studying the behavior of 
the total displacement Ar 12 rather than focusing only on 272, the component 
in the direction of Ar 01 . Ar 12 can be decomposed into the deterministic part 
((#12) given by Eq. 2 and (2/12) = 0), and a stochastic part. Both the determin- 
istic and stochastic parts may may depend on roi. To measure the importance 
of the stochastic part, we compute cry = (xf 2 ) — (xu) 2 and 07 = (yf 2 ) — (yu) 2 , 
shown in Fig. 4 by the connected symbols and unconnected symbols, respec- 
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Fig. 5. (a) The functions Pii(xi2|roi; t) and (b) -Pj_(yi2 l^oi? t), for <f> = 0.52 (a liquid), 
with t = t* = 600 s. The open circles are all data for roi < r cagc = 0.4 /im and 
the closed circles are for roi > r cagc . The gaussian fits are shown as dashed lines for 
7*01 < ^cage and solid lines for roi > r cagc , and have widths of a ps 0.13 /im for all 
except o"||(roi > r cagc = 0.22 fim. Similarly, the nongaussian parameter c*2 = 1.8 for 
all except -Py (xi2|roi > ?* cagc ), which has Q2 = 1.0. 

tively. The behaviors of the parallel and perpendicular components are similar 
at small values of r i, but differ markedly when the original displacement has 
a larger distance r m [11]. The transverse component <tj_ is nearly constant 
as a function of r i, but cry becomes much larger when r i is larger. Again, 
any dependence whatsoever on r i is indicative of memory in the system, and 
the increase in cry reflects a memory of mobility. Particles which move large 
distances originally (large values of roi) are more mobile subsequently (large 
values of cry), and in particular are more mobile along the direction of the 
original motion. 

Confirmation of this is seen by plotting the distribution functions -Py^^lroi) 
and P±_{yi2\ r m) m Fig- 5, where the open circles are for r i < r cagc and the 
closed circles are for r i > r cage . Gaussian fits to these distribution func- 
tions are shown by the lines. All of the functions appear similar, except for 
P\\(xi2',rm > r ca g e), which is significantly broader [solid circles in Fig. 5(a)]. 
Thus, particles which originally have larger displacements are more likely to 
continue moving in the same direction (large £01 > 0) or more likely to move 
a large distance backwards (x m < 0), but slightly less likely to stay in the 
same position. Moreover, as cry 7^ crj_, the particles undergoing cage rearrange- 
ments move in a highly anisotropic fashion. The distributions are broader than 
Gaussians, as can be seen by comparing the symbols to the lines; this is a re- 
flection of the underlying broad distributions of the displacements, as shown 
in Fig. 3(a). Note also that the distributions shown in Fig. 5 are symmetric 
about the peak; this is unsurprising for P^yu) and perhaps more surprising 
for P||(xi 2 ). 
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4 Discussion 



We have studied the microscopic motion of thousands of tracer particles in a 
concentrated colloidal sample, in order to understand the dramatic dynamical 
changes near the glass transition. In particular, near the glass transition, par- 
ticles are confined to transient cages, resulting in temporal anticorrelations in 
particle displacements. We find that caging can be described as a deterministic 
anticorrelated motion, plus a stochastic part. The deterministic part is due to 
memory provided by the caging particles, which must adjust their positions 
to allow a particle to move, and subsequently push that particle back toward 
its original position. By quantifying these effects (as given by Eq. 2), we can 
connect the properties of the cage directly to the subdiffusive growth of the 
mean square displacement (MSD), shown in Fig. 1(a). The connection is quite 
good, as seen by comparing the lines and symbols in Fig. 1(b). 

The long time behavior of the MSD is diffusive, as seen in Fig. 1(a). This can 
be thought of as due to the random walks taken by the individual particles, 
which alternate between being stuck in cages for a random duration, and a 
cage rearrangement motion of random length (see Fig. 2). A simple possibility 
which leads to diffusive motion at long times is that the cages responsible 
for the subdiffusive plateau have finite lifetimes with a characteristic time 
scale. An alternate possibility is that the cage rearrangement motions could 
be Levy flights. Levy flights are motions with an infinite mean square step 
size, in other words, cage rearrangements would involve movements that carry 
particles large distances. In such a way, diffusive motion at long times could 
be due to a competition between cages with infinite mean lifetime, and mo- 
tions with infinite mean square lengths [4]. (These possibilities would suggest 
that the distribution for cage times and/or step sizes are power laws, for ex- 
ample P(Ax) ~ (Ax)~ u for the cage rearrangement displacement Ax with 
1 < v < 3.) Levy flights seem possible when looking at the broad tails shown 
in Fig. 3(a). However, at best Fig. 3(a) shows a truncated Levy distribution. 
We do not see any particles making dramatic displacements much larger than 
their own radius; the trajectories shown in Fig. 2 making small adjustments 
are typical. It seems likelier that the characteristic step size is r cage , a small 
and finite distance, and thus the diffusive growth of the MSD as At — > oo is 
due to a finite cage lifetime [18]. In glassy samples, the cage rearrangements 
are no longer allowed, and thus the MSD will be subdiffusive at all times, and 
perhaps asymptotically reach a plateau; thus we expect these concepts to be 
even more useful in understanding the strange kinetics of nonergodic glassy 
samples. 

We thank B. Doliwa and H. G. E. Hentschel for helpful discussions, and thank 
A. Schofield for providing our colloidal samples. This work was supported by 
NSF (DMR-9971432) and NASA (NAG3-2284). 
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